Final Project

library(leaflet)
library (raster)
library(rgdal)
library(RColorBrewer)
library(devtools)
#devtools::install_github("rstudio/leaflet")
setwd("C:/Users/crish/OneDrive/Documents/Columbia SIPA/2017 Spring/Data Visualization/Final Project/")
library(readr)
iycf <- readstata13::read.dta13("data/MAP_2015_General_Info_IYCF.dta")
sam <- readstata13::read.dta13("data/2015_SAM_TOTAL.dta")
mnp <- readstata13::read.dta13("data/2015_MNP.dta")

Infant and Young Child Feeding (IYCF)

library(raster)
world_spdf <- shapefile("data/world_map/TM_WORLD_BORDERS_SIMPL-0.3.shp")
# Merge the QoG data to Shape Files
library(dplyr)
IYCF2 <- world_spdf@data %>% left_join(iycf, by = c(FIPS = "fips"))
world_spdf@data <- IYCF2

Make map…

leaflet(world_spdf) %>% 
    setView(lat=10, lng=0 , zoom=2) %>%
    addPolygons(
    stroke = FALSE, fillOpacity = 0.5, smoothFactor = 0.5, 
    color = ~colorNumeric("RdYlGn", Health_facilities)(Health_facilities),
    popup = paste("Country:",world_spdf$NAME,"<br/>",
                  "Health facilities:",round(world_spdf$Health_facilities,1),"","<br/>",
                  "Community health workers:",round(world_spdf$Community_health_workers,0),"")) %>%
  addLegend("bottomright", 
  pal = colorNumeric("RdYlGn", world_spdf$Health_facilities), values = ~Health_facilities,
  title = "Number of health care facilities", opacity = 0.5)

Micronutrient powders

Children given Micronutrient Powder, per region, 2015

childrenMNP <- c(203552,2285242,570648,4263234,446194,1785867,612016)
Region <- c("CEE/CIS", "East Asia and Pacific", "Eastern and Southern Africa", 
            "Latin America and the Caribbean", "Middle East and North Africa",
            "South Asia", "West and Central Africa")
MNP <- data.frame(childrenMNP,Region)
names(MNP) <- c("childrenMNP","Region")
library(ggplot2)
library(ggthemes)
library(scales)
ggplot(data=MNP, aes(x=Region, y=childrenMNP, fill=Region)) + geom_bar(stat = "identity") + 
  scale_fill_brewer(palette="YlGnBu") +
  ggtitle("Graph 1. Number of children given Micronutrient Powder, per region, 2015") + 
  labs(x="", y = "Children given MNP", size=8) +
  theme_classic(base_size = 12) + theme(text = element_text(color = "gray20"),
                                        axis.text.x = element_text(NULL),
          legend.position = c("bottom"), legend.direction = "horizontal",
          legend.justification = 0.05, legend.text = element_text(size = 9, color = "gray10"),
          legend.key.height=unit(1,"line"), legend.key.width=unit(1,"line"),
          axis.text = element_text(face = "italic"),
          axis.title.x = element_text(vjust = -1),
          axis.title.y = element_text(vjust = 2),
          axis.ticks.y = element_blank(),
          axis.line = element_line(color = "gray40", size = 0.5),
          axis.line.y = element_blank(),
          panel.grid.major = element_line(color = "gray50", size = 0.5),
          panel.grid.major.x = element_blank(),
          plot.margin = margin(t = 0, r = 0, b = 40, l = 5, unit = "pt"),
          plot.title = element_text(face = "bold", color = "black", size = 11, hjust=0.5))

#grid.text("Source: NUTRIDASH 2015", x = .01, y = .03, just = "left", draw = TRUE))

La misma pero interactiva…

library(plotly)
plot_ly(MNP, x = ~Region, y=~childrenMNP, color = ~Region, type = "bar", opacity=0.5)

plot_ly(iycf, x = ~Region, y=~Counselling_mothers, color = ~Region, type = "bar")

Ignoring 115 observationsIgnoring 115 observations

Scatter…

plot_ly(filter(edu_long, area_type=="county"), x = ~coll, y=~HS, 
                     color = ~year, type = "scatter",
                     text = ~paste("County: ", area_name, 'Year:', year)) 

Error in filter_(.data, .dots = lazyeval::lazy_dots(...)) : 
  object 'edu_long' not found

Severe Acute Malnutrition

library(readr)
iycf <- readstata13::read.dta13("data/MAP_2015_General_Info_IYCF.dta")